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We explore the morphological stability during the growth of strained multilayer structures in 
a dynamical model which describes the coupling of elastic fields, wetting effect, and deposition 
process. We quantitatively show the significant influence of the wetting effect on the stability 
properties, in particular for short-period multilayers. Our results are qualitatively similar to recent 
experimental observations in AlAs/lnAs/InP(001) system. We also give predictions for strain- 
balanced multilayers. 

Strained periodic multilayer films composed of different material layers have attracted much attention since such 
material structures can have tunable electronic properties. Among the actively investigated systems is the multi- 
layer with alternating tcnsile/compressive layers coherently grown on a substrate, resulting in a modulated structure 
like short-period superlattice [e.g., GaAs/InAsI or AlAs/In As^ii on InP(OOl)] or a multiple quantum well [e.g., 
GalnP/lnAsP on InP(OOl)] 5 . Generally morphological instability occurs in such coherent lattice-mismatch multilayer 
structures, driven by a gradual release of misfit-generated stresses when dislocations are absent. Consequently, undu- 
lations are observed in the strained layers, instead of ideally flat interfaces. This stress-driven morphological instability 
can lead to lateral compositional modulations in vertical short-period superlattices^Si&i as a result of layer thickness 
modulations (caused by interface rippling) and different material components in adjacent layers^^ The spontaneous 
formation of these lateral modulations is also a promising way to self-organize fabrication of low-dimensional quantum 
heterostructures, especially quantum wiresi 1 ' 2 ' 3 

Although understanding of the detailed growth mechanism in multilayer structures is important and of much interest 
for both fundamental studies and device applications, it is still far from complete due to the complexity of the system. 
This complexity arises from the coupling of strain fields in different layers, and the nonequilibrium nature of the 
growing film for which material deposition rates play an important role in the pattern formation. Previous theoretical 
analyses use some approximations to determine the elastic fields, treating each buried island^ or nonplanar interface^ 
of the multilayer as a misfitting inclusion in a semi-infinite homogeneous medium. Recently, a general method has 
been developed to directly calculate the elastic state of the multilayer system^ 

Wetting effects have been considered in the context of single-layer strained film growthj2iiS but not for multilayers. 
In this letter, we attempt to remedy this. The wetting effect arises from the change of material properties across 
interlayer interfaces within the system, and is especially important for thin composite layers. Here, we incorporate 
the wetting effect (arising from nonlinear elastic contributions 10 ) in the early morphological evolution of growing 
multilayer films. We focus on the stability properties of the system and on comparison with recent experiments. 

We consider a three-dimensional coherent multilayer structure with k stacking layers deposited on a semi-infinite 
substrate. In a general case, each composite layer i (1 < i < k) could consist of a material different from all the other 
layers and substrate, with average thickness U and height profile z = hi{x, y) = J2]=i h + S q ^t(l) exp(iq x x + i%y)- 
The film/substrate interface is described by /io(q), and the top strained layer k by a surface morphological perturbation 
fofe(q, t) and by an average thickness Ik = Vfet where Vk is the deposition rate. 

The basic state of the system is chosen to be a uniform film of flat growing surface and planar underlying interfaces, 
with all strains zero except for u zz .i = Ui = ej(l + Vi)/(1 — Vi) due to Poisson relaxation in the z direction. £j and 
Vi are the misfit and Poisson ratio of layer i, respectively. The elastic free energy of the multilayer is written as 
T c \ = d 3 r£Q + J2i=i Ih ' 1 d 3 r£i, where the clastic energy density for each layer i can be generally expressed as 

Si = £ i + S£i. Here £j is the energy density of the basic state and 5£i is the change in energy density due to a 
perturbation. Each of the two terms includes a linear contribution, obtained from linear elasticity, and a nonlinear 
contribution. The latter is often ignored, but is crucial for the wetting effect^ It leads to the dependence of £~i on 
layer thickness hi(x,y) — hi-\(x,y). This dependence is attributed to the coupling between the elastic relaxation of 
layer i and that of the underlying layer or substrate It was showniS that the thickness dependence of £i enhances 
the wetting of the substrate by the film in the case of a single layer, and consequently stabilizes the flat surface. 

For the energy density of the perturbed state, 8£i, we only keep the linear elasticity contribution 8£ \ m = cr ol p t i8u a p i i+ 
^a/3£,p8u a f3^Su^ p ^/2 (with o~ a p,i the stress tensor and \ a f3£ P the elastic modulus tensor), since the nonlinear contribution 
is of higher order in the perturbed strains Su a pj. Thus, we have £j = £ j + 8£ l ™. Using the results of Ref. for the 
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first-order surface elastic energy density, 5£ l k m , we obtain the first-order evolution equation for h k (q, t): 



dh k /dt = {v k -r kq 2 [d 2 f£l/dh 2 k ] hk=hk{t) }h k 
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with (7* = Tfe [i?'e^,(7 3 — 7fe(0)g 4 ] and / ol -(fti-i, fti) = i dz£i(hi — ftt-i, z). The derivatives of the latter quantity 
with respect to hi are taken with fixed hi—i. Here we assumed that surface diffusion is the dominant mass transport 
process, with negligible bulk and interlayer diffusion as well as frozen buried interfaces^ We used a step-function form 
for surface stiffness (7) and misfit stresses {cf X x,i — <?yy i) & t the interlayer interfaces^ In Eq. T k is the surface 
mobility and 7fc(0) is the surface stiffness at the orientation of the unperturbed surface. The average surface position 
is h k — Y^!i=i h + v kt, and E' = 2E(1 + v)/(l — v) with Young's modulus E and Poisson ratio v assumed to be 
identical for all the layers and substrate. Note that the corresponding results for single-layer film can be recovered 
from Eq. Q with k — 1. 

The evolution Eq. (JTJ is similar to that without the wetting effect^ showing explicitly the influence of interface 
morphologies and elastic properties of different buried layers on the top surface profile. The important results of this 
work arise from the term d 2 / e [ k /dh k , which represents the wetting effect. The exact solution of Eq. JTJ has a form 
similar to the solution in the absence of the wetting effect A The surface profile h k can be expressed in terms of the 
product of (2 x 2) characteristic matrices Lj (3 < j < k) as well as the profiles h\ and h 2 of first and second deposited 
layers. The elements of the bottom row of Lj are 1 and and those of the top row are Lj t i(q) and Lj^io), where 

Lj,i(q) = Hj(q,lj/vj) [l-T j E'q 3 e j (e :j - ej_i)Gj(g, Zj fa)] - L jj2 (q)/fJ,j-i(q,lj-i/v j -i), 

Lja{q) = -(J,j(q,lj/v j )e- ql i- 1 [r j e j G j (q,lj/v j )]/[r j -ie j - 1 G J - 1 (q,l j -i/v j - 1 )}. (2) 



Here 







In the following these general results, which can apply to both periodic and nonperiodic systems, are used to 
determine the morphological stability of periodic tensile/compressive multilayer structures, and in particular short- 
period superlatticesi2i2ii In this case the multilayer film is composed of alternating A and B layers, with their 
associated parameters such as nonzero misfits £* A i B y layer thicknesses l* A ^ B y surface mobilities ^* A i B y deposition 
rates v* A ^ and surface stiffnesses J A t B y Here we have expressed all the parameters in dimensionless form: l A — 
j A (0)/(E'e 2 A ) and v° A = T A E' 3 e 6 A /^ 2 A (Q), as well as e* B = e B /e A , T* B = T B /T A and 7* = j B (0)/j A (0). To evaluate the 
wetting effect, we assume the form df^Jdhj = £] m {I — Xj cxp[— — ftj-i)/ftML,i]} (for hj > ftj-i) appropriate for 
the single-layer case, (according to the qualitative modeling calculations in ref. IToL which yield this general functional 
form of wetting elastic energy for various lattice structures with diagonal bonds), where £] m = Ee 2 /(1 — v) is the 
linear elastic energy density of the basic state and /imlj is the 1 ML thickness of layer j (A or B). We also assume 
Xj = X — 0.05 for all strained layers. 

For this periodic A/B system, the characteristic matrices Lj presented in Eq. are reduced to two types, L^s), 
and the effective stability (for large number of grown layers k 3> 1) is determined by the maximal magnitude of the 
eigenvalues A of the matrix L^L^, as shown in Ref. 0. When max(|A|) > 1, the initial perturbations are amplified 
after further deposition, leading to a morphological instability; otherwise, the perturbations decay during the growth, 
leading to a stable flat film profile. 

Consequently, the stability diagrams for various growth and material parameters can be determined, as presented 
in Fig. ^ where in order to compare with experiments, the parameters are chosen to qualitatively represent those of 
the growth of AlAs/InAs/InP(001) short-period superlatticesi^ Note that the global strain of the whole multilayer 
film can be defined as cg = e A [l + (l B /l A )e B }/ (1 + Ib/Ia)- Thus for e* B = —I (as in Fig. 0, l B /l A — 1 corresponds to 
to = 0, while Ib/Ia > 1 or < I to tensile or compressive films. Similar to previous calculations without the wetting 
effect^ the stability diagram shows a clear asymmetry: the multilayer film is more stable when the global strain ec 
has the same sign as the misfit of the layers with smaller surface mobility. For the system of AlAs/InAs/InP(001), 
rinAs(Fyi) S> Taias^b), and then the globally tensile film {eq < 0) is more stable since eaiAs < 0, as seen in Fig. 
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^a) and in experiments. Furthermore, as shown in Fig. [Ja), without the wetting effect (dashed line), the flat film is 
stable only when the global strain is tensile. The situation is qualitatively different when the wetting effect is included 
[solid line in Fig. IHa)]: both the globally tensile and compressive flat films can be stabilized, although the asymmetry 
in the stability diagram is still obvious. 

The earlier result with the consideration of the wetting effect is qualitatively similar to experimental observations^ 
in the following way: We plot [Fig. Dfb)] max(|A|) as a function of global strain, eg, for v* A — 0.05, which corresponds 
to the horizontal dot-dashed line in Fig. ^a). To evaluate eq, we set ea = 3.2% according to the InAs/InP misfit. 
As explained earlier, when the maximal value of A is equal to 1, the multilayer film is stable, with a smooth surface; 
when it is larger than 1, the film is unstable, leading to a rough surface. Also, the larger the value of max(|A|), the 
rougher the surface that results, as obtained qualitatively from our theoretical results^ hk cc |A| fc//2 for k 3> 1. Thus, 
this figure can be compared with the recent rms roughness measurement of Norman et al.X where the surface is flat 
or with small roughness for samples under global tension or small compression, but shows a rapidly increasing degree 
of roughness for larger compression. This experimental finding can be qualitatively inferred from our result including 
wetting effect [solid line in Fig. ^b)], but not from the result without wetting (dashed line). 

We also find that with the wetting effect, the stability boundary value of I b /Ia, above which stabilization is possible, 
depends sensitively on l\ + l* B , but only very slightly on T* B . The stability boundary value is closer to 1 (the value in 
the absence of the wetting effect) when l* A + l* B is larger, showing that the wetting effect has less influence for thicker 
layer periods. 

Since some multilayer systems are grown under strain-balance conditions {eq — 0)^ it is interesting to study the 
corresponding stability properties. The diagrams of l* B /l\ (= — 1/e* b for eq — 0) versus v* A are shown in Fig. where 
the stability asymmetry remains and is independent of relative mobility T* B , indicating more stability for asymmetric 
structures^ From Fig. [21 one can see that with the strain-balance constraint, there is no qualitative change due to 
the wetting effect. Also, for short period (e.g., l* A + l* B = 0.3) and large mobility difference (e.g., T* B = 0.001), the 
wetting effect tends to render the film more stable. 

In summary, we have studied the interplay of elasticity, wetting layer effects and material deposition in coherent 
multilayer growth, and determined the stability and morphological evolution of growing films. The important influence 
of the wetting effect for short-period multilayers has been calculated, reproducing recent experimental measurements. 
We hope these results will be helpful for practical growth of low-dimensional nanostructures. 
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FIG. 1: (a) Stability diagram of 1*b/1*a vs rescaled deposition rate v\, with {\ = 0.05, solid line) or without (% = 0, dashed 
line) the wetting effect. The parameters are analogous to those of AlAs/InAs on InP(OOl) superlattices growth: T* B — 0.001, 
Ia + I'b — 0-3, es = — 1, Vg/v A — 1, and — 1. (b) The maximum value of |A| as a function of multilayer global strain eg, 
with the same parameters of (a) except for a fixed deposition rate v\ = 0.05, corresponding to the horizontal dot-dashed line 
in (a). 




FIG. 2: Stability diagram of 1*b/1*a vs v a m the strain-balanced condition, with (x = 0.05, solid line) or without (x = 0, 
dashed line) the wetting effect. The effects of different relative surface mobilities (F* B = 0.01 and 0.001) and multilayer periods 
[l A + l* B — 0.3 and 0.9 (inset)] are shown. Other parameters are v* B /v\ — 1 and j B — 1 as in Fig. 



